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We present a study of the high energy stability of a minimal complex singlet extension of the 
Standard Model with or without dark matter (CxSM). We start by obtaining the beta functions of 
the couplings of the theory from the effective potential and then perform the RGE evolution for the 
allowed parameter space of the model at the electroweak scale. We obtain the scale up to which 
the model survives and combine this information with all the LHC measurements as well as bounds 
from dark matter detection experiments as well as the relic density best measurement from Planck. 
We conclude that scenarios which solve both the dark matter and stability problems must contain a 
dark particle heavier than 50 GeV and a new visible state (mixing with the SM-like Higgs) heavier 
than 170 GeV. 


I. INTRODUCTION 

The recent success of the Higgs boson discovery at the CERN Large Hadron Collider (LHC) experiments [U 
[2], while putting the Standard Model (SM) on a firmer ground, is pushing us towards finding solutions to 
outstanding problems which are not yet addressed in the model. Such are the question of the nature of dark 
matter, what is the mechanism generating the observed matter anti-matter asymmetry or even what solves the 
apparent metastability of the SM Higgs potential [3] . 

A partial answer to these questions may be about to appear in the next runs of the LHC, and could manifest 
itself in a minimal form. An example of a minimal parametrisation which naturally connects to the Higgs sector 
under scrutiny at the LHC is given by the scalar singlets framework. This provide, in addition, quite a natural 
connection into hidden sectors [4]. In fact, in the SM, there are only three gauge singlets that we can build with 
dimension less than four. In particular, the only scalar is H {H denotes the Higgs doublet), which means 
that scalar singlet fields naturally couple to the Higgs through this operator 

= Vsm{Wh) + Wh X + VneMi) • ( 1 ) 

Here we split the potential in three parts: i) the SM-like part, Vsm; H) The SM Higgs singlet operator coupling 
to some operator, with mass dimension up to 2, which depends on the new singlet fields (j)i, and iii) the 

new (purely) singlet potential L^ew this class of models, the physical mass states couple to the SM only 
through the Higgs field fiuctuation, if they mix, otherwise they are dark particles. Their mixing factors (/^a, 
a = 1,..., #scalars) suppress their couplings to other SM particles compared to a SM like Higgs coupling. At 
tree level, if we denote the SM Higgs boson fiuctuation by h and the mass eigenstates by Ha, then we have a 
sum rule for the Ka [5] 

= , ^|«aP = l. (2) 

a a 

In this study we focus on complex singlet (i.e. two scalar) singlet extension which provides both a simple 
parametrisation of dark matter (through the state that do not mix with the SM Higgs fiuctuation) and can be 
made compatible with electroweak baryogenesis as to generate the matter anti-matter asymmetry IHH]. We 
will find that, in addition, such model can also improve the stability of the SM at high energies m- 
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II. A MINIMAL COMPLEX SINGLET MODEL WITH DARK MATTER 

We consider a model where the SM is extended by a complex scalar singlet § = (5'+iA) /^/2 with potential [11] 

1 / = |§|2 + || g |2 ^ ^| g |4 ^ ^ ^ ^ _ ( 3 ) 

This can be seen as a softly broken version of a U(l) symmetric model which preserve a Z 2 parity for the 
imaginary part, A. This minimal complex singlet model contains two possible phases which break electroweak 
symmetry as necessary in the SM. If we define the vacuum expectation values, Vi, of the new fields through 
S = Vs -\- s, A = va P while keeping the SM Higgs vev v then we have: 

• Z 2 phase {vs 7 ^ 0, = 0) - symmetric phase: Here the two fields h, S mix to form the two visible mass 

eigenstates and A is a dark matter candidate. We highlight ki 26 = cos a and /^new = sin a in the mixing 
matrix 

( ^126 \ / 1^126 -sina 0\ f h\ 

-Hnew 1 ~ I ^new COS CX 0 1 I 5 1 . (4) 

Hbm ) \ 0 0 ^/V^/ 

• JPf phase {vs 7 ^ 0, 7 ^ 0) - spontaneously broken phase: All scalar massive degrees of freedom are visible 

and mix. Again we highlight the Ki in the mixing elements Rij 

( Hi 26 \ I 1^126 Rls Rla 

-Hlight I ~ I flight R 2 S -^ 2 a 
-Hneavy J y ^Heavy -^3s -^3a 

In both phases we match the first visible scalar to have the mass of the recently found Higgs. In the broken 
phase, the other two are ordered in mass relative to each other (but they can both be lighter or heavier than 

the Higgs or one lighter one heavier than the Higgs). This minimal model is representative of all kinematically 

interesting situations in terms of visible and invisible channels in a sector with two new real singlet degrees of 
freedom at the LHC. 

Finally it is important to note that, in this model, observables such as for example cross-sections, can be 
easily related to SM observables through the Hi factors (for example the production cross section for one of 
these scalars is cr{Hi) = nf(TsM{Hi) - see also (Hj). 



A. Renormalisation group equations and stability 


In this study, our perspective was to perform a global scan over all possible scenarios combining the phe¬ 
nomenological constraints, at the electroweak scale, with the requirement that the theory remains stable up to 
a high energy scale (such as the GUT or Planck scale). Previous global scans for this model where mostly tree 
level mm, and one-loop studies focused on particular points/scenarios [13]. 

To obtain the evolution equations for the couplings of the theory, to run their input values from the electroweak 
scale to high energies, one needs to compute loop corrections of the following two quantities (to keep track of 
the loop order, we define 5 = fi/(167r^)): 


^ 14^^ = y(0) + -h 52 u(2) _ 


^=0 




Gp = uy^+eni;> + e^u'y> + .. 


r(i) 


^ 2 tt ( 2 ) 


divs. 


dt 

1 dvj 
. Vi dt 


+ ... 


+ • • • 


(6) 


In the remainder t = log /i, where pt is the renormalisation energy scale. On the first line we indicate the loop 
expansion of the effective potential, iAff, with the n-loop order correction. The scale invariance of the 
effective potential provides the beta functions which determine the evolution equation of each coupling (all 
couplings denoted collectively by L here). Using the general formalism in [14] we have computed them at two 
loops. 

The second quantity to compute (second line of Eq. (HI) is the inverse propagator (or two point function) 
which encodes the anomalous dimensions ) • These are responsible for the evolution of vacuum expectation 
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values of fields (right hand side). Again, this is given by the tree level inverse propagator plus loop corrections, 
the self energies and the anomalous dimensions are extracted from their divergences m- 

In summary we obtained the Renormalisation Group Equations (RGEs) at two loops including all SM con¬ 
tributions m- 

To provide low energy input for the RGE running let us consider the loop expansion of couplings (given 
physical input such as masses and VEVs) and corresponding beta functions at the Z-boson mass scale: 


L{M,) = +sW'>imlvi) + ... 

/?L = + ... 


( 7 ) 

(8) 


It is clear from these expressions, that if we trunc ate the coupling expansion at tree level (one-loop) the beta 
function is of one-loop (two-loop) order. In Sect. Ill we evolve the couplings to high energies using one and 
two loop RGEs which we have checked with tree level and one-loop accurate low scale input data [ 12 ] . In such 
study, in addition to imposing a correct electroweak symmetry breaking pattern and stability at the low scale, 
we also impose that no runaway directions are developed and that perturbative unitarity holds, at all scales. 
This second condition, in practice, also prevents Landau poles. Note that strict stability up to the high scale, 
which is what we are requiring, is not possible in the SM as shown in [3]. 


B. Phenomenological constraints 

In all scans we have used the Scanners code [T6|, where a model class was implemented to impose all 
theoretical constraints (vacuum stability, boundedness from below and perturbative unitarity). In the phe¬ 
nomenological scans we have implemented constraints from electroweak precision variables (STU) directly in 
Scanners, and collider constraints were applied using HiggsBounds (to set 95% G.L. on new unobserved scalars) 
and HiggsSignals to obtain the probability of the model point to fit the data [13 [18]. All SM-like decay widths 
were computed with Hdecay [19]. Dark matter constraints from the Planck data [20] for the relic density and 
limits on direct detection cross section from LUX [21] were applied by computing the corresponding quantities 
for this model with micrOMEGAS [22] (see also [12] for details). 


III. DISCUSSION OF THE RESULTS 

In Eig.|III[ we show allowed parameter space points where only theoretical constraints were imposed (vacuum 
stability, boundedness from below and perturbative unitarity). The colour code of all panels corresponds to the 
energy scale at which either a runaway direction develops or perturbative unitarity is violated. On the top left 
panel, in the vertical axis we have the Higgs quartic coupling versus the mass of the new visible scalar mixing 
with the Higgs, whereas the right panel contains the coupling suppression of the Higgs found at the LHG. The 
two boundaries on the left panel are due to the minimum conditions and we can see that there is a yellow band 
for 0.5 < A < 1 for which the model is stable up to the Planck scale which we call the stability band. On the 
right panel we find that if we impose stability up to an intermediate scale, 10 ^^ GeV, there is a lower bound on 
the new mixing scalar mass of roughly 140 GeV. 

In the bottom panels we present, for the broken phase, the same coupling versus the mass of the lightest (left) 
and the heaviest (right) of the new visible Higgses. The two boundaries of the stability band of the dark phase 
are now mapped to the upper and lower boundaries of the stability bands of the lightest and the heaviest of the 
new scalars, respectively. What is interesting to note, is that: i) if a state lighter than 140 GeV is found, the 
dark matter phase is disfavoured and ii) again, we need one new scalar mixing with the Higgs which is heavier 
than 140 GeV to take the SM from a metastability scenario into stability. 

Einally, in Eig. |HI[ left panel, we present the effect of the phenomenological constraints alone (without RGE 
running) on the (/^i 265 ^i^new) for the dark phase. All points are consistent with the Higgs signals at the 

LHG within 3cr. In the coloured points we have imposed that the new visible scalar cannot decay invisibly to 
the dark matter candidate (in contrast with the gray points which do not have this cut). In the right panel, we 
display a projection with the two new scalar masses (dark matter, tudm and where we impose stability 

up to the GUT scale, and that the model saturates the dark matter relic density measured by Planck (together 
with the limits from LUX). This pushes the 140 GeV lower bound on the new heavy scalar to roughly 170 
GeV [24] whereas tudm ^ 50 GeV. 
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FIG. 1: Parameter space projections displaying the with the scale up to which the theory remains stable (in colour) for 
the dark phase (top panels) and the broken phase (bottom panels). (Adapted from |12]b 
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FIG. 2: Left: SM-like Higgs coupling versus the mass of the new visible particle (in the dark phase) with phenomeno¬ 
logical constraints at the weak scale applied, consistent within la, 2a and 3a of the LHG Higgs measurements. Right: 
Gombination of all phenomenological constraints with the requirement of stability up to the GUT scale and the correct 
relic density. (Adapted from [l2]). 


IV. CONCLUSIONS 

We have obtained the two loop RGBs for a complex singlet extension of the SM with dark matter, which 
may help us improve our understanding of baryogenesis and the stability of the SM at higher energies, and that 
contains a rich LHC phenomenology (in terms of visible and invisible decays of two new scalar singlet degrees 
of freedom). In our stability analysis we found that both phases of the model have a stability band where the 
theory is better behaved at high energies. We also found that, for this to occur, the scalar spectrum must 
contain a new heavy scalar mixing with the Higgs with mass larger than 140 GeV. 

When combined with phenomenology, imposing that the model explains dark matter fully and that it stabilises 
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the theory, the lower bound on the new heavy scalar is pushed up to 170 GeV and the dark particle must be 
heavier than 50 GeV. 
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